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The diffusion Monte Carlo technique is used to calculate and analyze the excitation spectrum of 
a single 3 He atom bound to a cluster with N 4 He atoms, with the aim of establishing the most 
adequate filling ordering of single-fermion orbits to the mixed clusters with a large number of 3 He 
atoms. The resulting ordering looks like the rotational spectrum of a diatomic molecule, being 
classified only by the angular momentum of the level, although vibrational-like excitations appear 
at higher energies for sufficiently large TV. 
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I. INTRODUCTION 

The study of isotopic 4 He- 3 He mixed clusters is at- 
tracting a growing interest in recent years. From the 
experimental viewpoint, the diffraction of clusters from 
a transmission grating |]J has opened new perspectives in 
the detection and identification of small clusters. There 
is no evidence for the existence of the dimcr 4 Hc 3 Hc, 
but clusters 4 He/v 3 He with N = 2,3,4 have been defi- 
nitely detected 0] . It seems possible at present to resolve 
clearly the clusters of mass up to about 25 amu 0. As 
the weakly van der Waals He-He interaction is isotope 
independent, the properties of such mixed clusters are 
determined solely by quantal effects, namely the differ- 
ent zero-point motion and the different statistics of the 
two isotopes. It turns out that helium clusters are weakly 
bound systems, and the lighter ones are challenging for 
microscopic theoretical methods. 

The stability of small mixed clusters has been the ob- 
ject of several recent microscopic studies. Guardiola and 
Navarro have investigated clusters containing up to eight 
4 He atoms and up to 20 3 He atoms, based on both a 
Variational Monte Carlo (VMC) wave function PjL Q 
and the diffusion Monte Carlo (DMC) method |(| in 
the fixed-node approximation. Bressanini and collabo- 
rators 0, @ have considered clusters with up to 17 
4 He and up to three 3 He atoms, by means of the DMC 
method, also in the fixed node approximation. 

The DMC description is based in an importance sam- 
pling wave function which plays a triple role: it controls 
the variance of the ground-state energy, it carries on the 
quantum numbers and other properties of the consid- 
ered cluster, and it specifies the nodal (or set of nodal) 
surfaces. In particular, in Refs. U- lj L| the antisym- 
metry required for fermions has been taken into account 
by means of two Slater determinants, one for each spin 
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orientation, which have been built up in terms of har- 
monic polynomials in the Cartesian coordinates of each 
fcrmion. Moreover, a harmonic-oscillator (HO) ordering 
of the fermionic shells has been assumed. Although this 
seems a reasonable hypothesis, supported by the find- 
ings of density-functional calculations in medium sized 
mixed droplets |l0j . from a microscopic point of view 
there are no conclusive a priori arguments in favor of 
such an ordering. For instance, to describe mixed systems 
with three 3 He atoms, a configuration with total angu- 
lar momentum L = has been assumed in Refs. 
which corresponds to the filling of the single-particle lev- 
els ls 2 2s. In contrast, the \s 2 \p HO ordering, which has 
been assumed in Refs. 0,|U,Eli results in an angular mo- 
mentum L = 1. The comparison of respective binding 
energies indicates that the L = 1 state has a lower en- 
ergy than the L — one. In conclusion, a more general 
criterion to select the shell ordering needs to be specified. 

The aim of this paper is to determine the excitation 
spectrum of a single 3 He atom bound to a 4 HeAr clus- 
ter. The ordering of these single-particle levels will be 
relevant to describe mixed clusters with a higher number 
of 3 He atoms. It is also worth noticing that the knowl- 
edge of the one-fermion spectra in terms of the number 
of bosons is also relevant to determine the constants en- 
tering the rate equations which establish the formation 
chemical process [llj, an d thus the abundances in the 
production experiments. Our calculations are based in 
the DMC method, using an importance sampling wave 
function which carries out the orbital angular momentum 
of the 3 He relative to the 4 HeAr cluster. The DMC proce- 
dure is thus adapted so as to determine the lowest energy 
state of the subspace of orbital angular momentum L. In 
order to obtain the excited states within each subspace 
of angular momentum L we use an optimized form of the 
upper bounds provided by the sum rules method. 

The paper is organized as follows. In Section II we 
briefly review previous investigations of 4 Hejv 3 He clus- 
ters. In Section III we present a detailed description of 
the method used to study the ground and excited states 
of these clusters. Our results are presented and discussed 
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in Section IV. Some final comments are given in Section 
V. 



II. A SURVEY OF PREVIOUS RESULTS ON 



atoms, so that the cluster 4 Heiv 3 He can be considered as 



HEat^HE clusters 



A pure 4 Hejv+i cluster is described by the Hamiltonian 



H 



2 N+l 



2m 



(1) 



where is the mass of 4 He. and V(r) is the interaction 
potential. Recent forms [la . Ibll ITq of this interaction 
are quite similar and we will use along the paper the one 
known as HFD-B potential. 

Given that the He-He interaction is a consequence of 
the interaction between the electrons in the atoms, it is 
independent of the mass or the spin of the nucleus. To 
convert the (iV + l)-th atom into an 3 He atom, thus deal- 
ing with the cluster 4 He,/v 3 He, corresponds to a simple 
change in the Hamiltonian 



H 



ti 2 

2TO4 



N+l 



E^ 2 
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N+l 



E v ( r a)> ( 2 ) 
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where 7713 is the mass of 3 He. The corresponding many- 
body problem is then not much different from that of a 
pure 4 He cluster. In the rest of the paper we will use the 
subindex F instead of N + 1, to alleviate the notation. 

Mixed 4 Hejv 3 He clusters containing a single fermion 
have been investigated using several theoretical methods. 
The hrst systematic study of the excitation spectrum of 
the 3 He atom was made by Dalfovo ^(|, based on a zero- 
range density functional. The use of a non-local finite- 
range density functional 01 results in small quantitative 
differences, related to the fact that finite-range function- 
als are more repulsive than zero-range ones. As the size 
on the drop increases, the 3 He atom is pushed to the sur- 
face region, due to its large zero-point motion, and for 
large enough clusters, the centrifugal term L(L + l)/r 2 
entering the Schrodinger equation for the 3 He atom can 
be treated as a perturbation. Actually [l^. the gen- 
eral trend of the spectrum is a series of tight rotational 
bands on top of radial excitations, related to the num- 
ber of nodes in the radial wave function. In the limit of 
very large N, the spectrum of the 3 He atom becomes in- 
dependent of the quantum angular momentum, forming 
the analogous of the two-dimensional Andreev states in 
bulk helium. 

The so-called Lekner approximation was used in 
Ref. 0], where a VMC calculation was performed, as 
well as in Ref. based on the Hypernetted-Chain 

method in its optimized version. Such approximation, 
used by Lekner [ljj to analyze the Andreev states in bulk 
liquid, assumes that the pair correlations between 3 He 
and 4 He atoms are the same as those between pairs of 4 He 
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a perturbation of the cluster HeAr + i, the perturbation 
being given by 



H, = 



( & 



\2r7l4 2?Ti3 



(3) 



Further elaboration of the perturbation scheme results in 
a single-particle Schrodinger equation describing the 3 He 
atom with an effective potential given by 



h 2 V 2 7MQ 
2m 3 y/p 4 (r) 



(4) 



where pa{t) and T±{r) are respectively the 4 He particle 
and kinetic energy densities, both defined in the unper- 
turbed system with N + 1 bosons. Notice that the Lapla- 
cian operator acting on the square root of pn produces 
a strongly attractive force peaked at the surface of the 
cluster. 

In order to classify the resulting spectra, Krotschcck 
and Zillich Il8l defined an effective wave number k = 



\/L(L + l)/R for each excitation characterized by the 
orbital angular momentum L, where R is the equivalent 
hard sphere radius R = ^/5/ 3?' rm s , defined in terms of 
the root mean square radius r rms of the droplet. By plot- 
ting all excitation energies as a function of k for a large 
number of clusters, Krotscheck and Zillich found that all 
results fall reasonably well on a universal quadratic line, 
in nice agreement with the density-functional results. 

Both sets of results may be approximately pictured 
as molecular rigid rotor, in which the 4 Hejv 3 He cluster 
is viewed as a two body system, the cluster formed by 
the N bosons plus the single fermion, tied by a spring 
with a rather large rigidity constant. Various angular 
momentum L states are associated to each vibrational- 
like state of the spring, obeying the law 



5E L « h 2 L(L+ 1)/2J = KL(L+ 1). 



(5) 



This equation defines the rotational constant K (units of 
energy) in terms of the momentum of inertia I, the later 
being proportional to R 2 , where R is the average dis- 
tance of the fermion to the center-of-mass of the bosonic 
cluster. 

In the case of light clusters (N < 40) the calculations 
of Ref. [Tsj find appreciable deviations from the univer- 
sal behavior Eq. JSJ. Small 4 HeAr 3 He clusters have been 
studied in Ref. pfl based on the DMC method, and in 
Refs. iM without the use of the Lekner approxima- 
tion. To this respect it is worth stressing that the Lekner 
approximation is basically a weak coupling description of 
the interaction of a 3 He atom with a 4 He cluster, because 
it does not include the perturbation that the outer 3 He 
atom should generate on the binding cluster. One may 
expect this picture to be satisfactory for large bosonic 
clusters, but inappropriate for small N. As a conse- 
quence, the interesting area to explore corresponds to 
the case of a small 4 He cluster, which may be apprecia- 
bly modified by the 3 He atom. 
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III. GROUND AND EXCITED STATES 
DESCRIPTION 

To carry out the Diffusion Monte Carlo calculation re- 
quires an importance sampling or guiding wave function, 
which incorporates as much as possible the characteris- 
tics of the system to be described. In particular, due 
to the strong short-range repulsion of the atom-atom in- 
teraction, it is advisable to introduce at least two-body 
Jastrow correlations. Moreover, the guiding wave func- 
tion must confine the system and, finally, it has to include 
the bosonic symmetry related to the 4 He atoms and the 
desired quantum numbers for the 3 He atom. In this Sec- 
tion we shall describe a variational wave function which 
will be used to calculate the lowest-energy states for a 
given value of the angular momentum L of the system. 
Excited states corresponding to radial excitations will be 
estimated by means of sume-rule techniques. 

A. Ground state and angular momentum 
excitations 

In order to describe the system 4 HeAr 3 He in a state 
where the 3 He atom is in an orbital angular momentum 
L with respect to the 4 HeAr system, a simple but never- 
theless complete wave function is given by 

tt(ri,...,rjv,r F ) = $ s (ri, . . . , rjv) (6) 
$M(ri, . . . ,rjv,r F )$ L (r F - R B ) , 

where the subindexes B and F stand for bosons and 
fermions, respectively, whereas M refers to the mixed 
boson-fermion part of the wave function. The bosonic 
coordinates run from 1 to TV and the coordinate of the 
fermion is labelled by F. Finally, R# represents the 
center-of-mass coordinate of the bosonic subsystem. This 
model wave function includes an internal bosonic part 
(&b) and the coupling of the fermion to individual bosons 
(3>m) as we H as t° the bosonic cluster 
We have taken and <I>m to be of the Jastrow form 

N 

^> B (r 1 ...r N ) = J] e^> (7) 

i<j=l 
N 

^(ri.-.r^rf) = ]Je^ riF \ (8) 

i=i 

with 

/s.M(r) = ~ ~Pb,m t - ( 9 ) 

The two-body correlation terms include a short-range 
part, associated with the parameters bs and 6m, and 
a long-range confining part associated with the parame- 
ters ps and pm- The short-range part is mainly related 
to the small-distance behavior of the relative two-body 



wave function. Consequently, the parameters 6_b,m and 
v have been taken to be the same for all systems stud- 
ied, and in our calculations they have been kept fixed to 
the values bs = 2.95 A, 6m = 2.90 A and v = 5.2, as 
obtained in our prev ious calculations for pure 4 He and 
3 He clusters pfl l2lL l22l |. by direct minimization of the 
expectation value of the energy. On the other hand, the 
long-range confining parameters pb,m have been deter- 
mined by means of the VMC method. 

In the absence of the last term of the importance 
sampling wave function © we would describe a state of 
null angular momentum, explicitly translational invari- 
ant and including the bosonic symmetry of the indistin- 
guishable 4 He atoms. The role of the last term of 
Eq. ©, describing the motion of the fermion, is to de- 
termine the value of the orbital angular momentum L. 
It has been taken as a long-range wave function depend- 
ing on the relative coordinate of the fermion with respect 
to the center-of-mass of the bosons, r = r F - R^. The 
explicit form used is the harmonic polynomial 

$ L (r) =r L P L (cos 9 r ). (10) 

This function is particularly simple and corresponds to 
a state with orbital angular momentum L and null third 
component. One could have taken a more sophisticated 
form, by putting a radial dependence different from the 
simple r L , but it is reasonable to expect the DMC al- 
gorithm to be able to improve this simple and compu- 
tationally convenient form. By making this part of the 
trial wave function to depend on the relative distance 
of the fermion to the center-of-mass of the bosons the 
translational invariance of the importance sampling wave 
function is not spoiled. Notice that if we had considered 
the function to depend on the distance of the fermion to 
the center-of-mass of the full system, the only difference 
would have been a trivial scale factor. 

With the structure of the importance sampling wave 
functions one may describe the lowest-energy states for 
each angular momentum L. It should be mentioned that 
apart from the case L — for which = 1, all other 
cases correspond to functions with a nodal surface, with 
nodes depending only on the angular variables. This fact 
has to be taken into account when using the DMC algo- 
rithm. 



B. Radial excitations 

When considering a subspace of angular momentum L 
the DMC algorithm gives only the energy of the ground 
state of that subspace, and there is no information about 
the excited states of the same angular momentum. A way 
to have an estimate, actually an upper bound, of the first 
excited state is to use the sum rules method [23| . 

Consider the exact ground state for a given angular 
momentum L, represented here by ^ol, and the full set of 
eigenstates of this subspace ordered by increasing energy 
and represented by {f„i,B n; i}, n = 0, 1 . . .. Let Q(R) 
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be an arbitrary Hermitian operator which may depend 
on all atomic coordinates, which is assumed to be scalar 
under rotations, i.e., to commute with L and S. Let us 
consider the sum rule of order p 



M^[Q] 



£ol) p I0MQI*ol}| 2 . 



(11) 



where the sum extends to all eigenstates of the Hamilto- 
nian but the lowest energy state of angular momentum 
L. This is important because in order to obtain easily 
computable properties it will convenient to use the com- 
pleteness relation. Because of the assumed properties of 
Q only states with angular momentum I = L will con- 
tribute to the sum. The p = 1 rule fulfills the property 

m£\q] = J2( E nL - ^ol)|(*„l|Q|*ol)| 2 

> (E 1L -E 0L )M^[Q], (12) 

from which one obtains an upper bound to the energy of 
the first excited state of the subspace L 



A/f[Q] 



(13) 



The evaluation of the sum rules is simpler than seem, 
because of the relations 

M< 0) [Q] = (v& 0L |Q 2 |* 0L )-|(vI/ 0L |Q|* 0L )| 2 (14) 



M { L 1] [Q] = -{* 0L \[Q,[H,Q]}\V OL ) 



(15) 



The double commutator may be simplified for a general 
Hamiltonian of the form given in Eq. ©, obtaining 



Q. Unrestricted minimization will give rise to the unprac- 
tical relation Q\^ol) = I^il); its solution being equiva- 
lent to the solution of the many-body Schrodinger equa- 
tion for the excited state. An alternative is to optimize 
the operator inside a restricted subspace, which is the 
approach followed by Chin and Krotscheck 0, H3 and 
is closely related to the procedure of Krisna and Wha- 
ley 0. 

Here we have followed a simpler procedure, based in 
the linear expansion of the operator in a basis of easily 
computable operators. To determine the basis we have 
assumed a single-particle-like form for the operator, by 
considering that it depends only on r F — Rb, i.e., on the 
coordinate of the fermion referred to the bosonic center- 
of-mass 



Q = Q{r F -R B ). 



(17) 



This simple form preserves the translational invariance 
and does not spoil the boson symmetry of the 4 He sub- 
system. For this type of general operator Eq. (|16fl can 
be further simplified, since the following relation 



V,Q(r F -R B ) 



N 



V F Q(r F - R B ) 



holds. The resulting sum rule Mi becomes 



Mi = — <*oi||V F Q| 2 |tf i>, 

2jJL 



(18) 



where fi = (N m^rns) / (Nm^ + 7713) is the reduced mass 
of the 3 He atom and the 4 He cluster. 

In the calculations to be described below the monopole 
operator has been optimized by using a simple functional 
form depending on few parameters, 



M« = -2- 

L 2m 4 



JV 



(16) 



Q(r F - R s ) = E C m1m{^F - Rb) 



(19) 



m— 1 



2m 3 



(4-ol||V f Q| 2 |* ol ) 



with 



Note that to compute these expressions one only requires 
knowledge of the ground state wave function of the an- 
gular momentum L subspace. 

This method was used in Ref. |2^| to obtain upper 
bounds to the first L = excitation, as well as to the 
low-lying even-L states. Given that we are obtaining the 
L^O excitations directly from the DMC procedure, the 
sum rules method will be used here to obtain the ener- 
gies of the first excited states in each L-subspace. In the 
Appendix we use the sum rule method to also estimate 
L ^ excitations based in the knowledge of ^00 by re- 
laxing the scalar character of Q as an alternative to the 
direct DMC calculations. 

The upper bound given by Eq. I|13|l is a functional of 
the operator Q, so it may be variationally optimized by 
equating to zero its functional derivative with respect to 



g TO (R) = |r F - Ri 



(20) 



The minimization of the upper bound of Eq. (|13l) with 
respect to C m , for angular momentum L, gives rise to a 
generalized eigenvalue problem, 



with a Hamiltonian-like matrix 



(21) 



■ M mn = ^ (*0i|V F <?mV F <?n|*0L) (22) 

and a normalization matrix 

M° mn = (*oL|? m « B |*o£) (23) 

- (^0L\q,n\^0L)(^QL\qn\^QL)- 
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The lowest eigenvalue of Eq. (|21() provides an optimized 
upper bound. By inserting Eq. I|20[l into Eqs. <|22|) and 
(I23JI the matrices are further simplified to 



2fi 



mn(^ 0L \ \r F - R B 



m+n— 2 



I*ol) (24) 



and 



M<® = (*o L ||r F -R B | m +"|* OL ) (25) 
- (*ox| \tf - RbH*ol><*ol| |r F - R B | n |*oL>- 

It is worth keeping in mind that in our DMC calculations 

(p) 

the matrix elements Mmn are based in a mixed estimate, 
so the strict variational character of Eq. (|13|1 may be lost. 

Some insight on the structure of these excitation may 
be drawn by considering the leading m = n = 1 terms 
of both the Hamiltonian and norm matrices. For this 
one-dimensional subspace the upper bound to the radial 
excitation energy is given by 



Eil — E l 

h 2 



< 



(26) 



2 M (0L\ \r F - R B | 2 |0L) - |(0L| |r F - R B | |0L)| 2 ' 

The size of this bound depends on the difference be- 
tween the mean square radius and the squared mean ra- 
dius of the fermion with respect to the boson center-of- 
mass. This difference will be small if the distribution of 
the fermion with respect to the boson center-of-mass is 
sharply peaked near a given value say i?o, corresponding 
to a rather rigid spring. Then the denominator will be 
small and the radial excitation will have large energy. 

The excited states considered in this section result from 
excitations related to the distance between the fermion 
and the center-of-mass of the cluster of bosons. The un- 
derlying optimal wave function (not determined, how- 
ever, from the DMC procedure) should have at least a 
node along this radial coordinate, in order to be orthog- 
onal to the angular momentum L ground state, and may 
be properly termed as a radial excitation. 



C. Computational details 

The DMC algorithm [H ^ is nowadays a well- 
known and used technology. It is based in integrating 
the imaginary-time Schrodinger equation for an auxiliary 
function /(R, t) = $ var (R)5'(R, t) which is the product 
of a trial wave function 4> var and the true ground-state 
wave function ^(R, t). The solution is given in terms of 
an approximate small-time Green funtion G(R, R', r) 



/(R',t + r)= J dRG(R',R,r)/(R,<) 



(27) 



by means of a series of small time steps r. We have used 
the 0(t 3 ) approximate Green function [2||[3(| which pro- 
vides an 0(t 2 ) approximation for the energy. In our 



calculations we have used a set of 1000 walkers, on the 
average, a value of r = 0.0002 K -1 , 8000 iterations to 
settle down the system and 80000 iterations to compute 
the averages. For N > 40 we have used a smaller time 
step r = 0.00015 K -1 and the number of iterations has 
been doubled. 

The wave function is not definite positive when L =/= 0, 
and this leads to the known and irritating sign problem 
of the DMC algorithm. We have used the so-called fixed 
node approximation. The auxiliary function /(R, t) will 
remain positive if both functions $ var and W have (at any 
time) the same nodal surfaces. The fixed node approxi- 
mation consists in killing any walker which att emp ts to 
cross a nodal surface. It has been shown [23, |28l ISJ that 
this procedure leads to an upper bound to the ground 
state. 



IV. 



ENERGETICS OF 4 HEjv 3 HE CLUSTERS 



We have first performed a VMC calculation to deter- 
mine the free parameters pb and pm of the importance 
sampling wave function. As mentioned in the previous 
Section, the parameters & B , &m and v are fixed by the 
atom-atom interaction at very short distances, and do 
not depend on the size of the cluster. On the contrary, 
the long-range parameters pb,m are very sensitive to the 
number of bosons. Their values, determined by minimiz- 
ing the ground state expectation value of the Hamilto- 
nian, are reported in Tabled 



TABLE I: Optimal values of the parameters (in A -1 ) for the 
variational description of the system 3 He 4 Hejv, as a function 
of the number of bosons N. The last two columns contain 
the ground state VMC and DMC energies (in K). Figures in 
parenthesis are the standard deviations of the Monte Carlo 
calculations. 



N 



PB 



Pm 



Evmc 



Ey)mc 



5 0.1250 0.0900 
10 0.0672 0.0395 
15 0.0491 0.0321 
20 0.0365 0.0226 
25 0.0298 0.0187 
30 0.0253 0.0160 
35 0.0220 0.0140 
40 0.0195 0.0125 
45 0.0176 0.0113 
50 0.0160 0.0104 



-1.59(2) 
-8.99(6) 
-19.17(13) 
-31.09(12) 
-44.51(16) 
-59.2(2) 
-74.5(3) 
-90.5(3) 
-106.1(3) 
-122.8(5) 



-1.862(4 
-9.763(9 
-21.39(2 
-35.42(3 
-51.20(3 
-68.03(5 
-86.13(6 
-104.90(9 
-124.63(5 
-144.67(5 



One can see that the values of these parameters de- 
crease when the number of bosons increase, as it should 
correspond to a drop, its size growing with the number 
of constituents. Notice that, for any given number of 
bosons, the parameter pm , controlling the boson-fermion 
distance, is significantly smaller than the parameter p B 
controlling the boson-boson distance. This reflects the 
fact that the light particle stands near the surface of the 
bosonic drop. 
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FIG. 1: Raw DMC results (boxes with error bars) and the 
least-squares fits. From top to bottom, the figures show the 
dissociation energy and the excitation energies for L = 1 — 4, 
all measured in K, as a function of the number of bosons 
N. The dashed line in the lower four figures is the smoothed 
dissociation limit. 
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In Table [I] the ground state energies obtained within 



the VMC calculation as well as with the improved DMC 
method are also displayed. The later are significantly 
lower than the former, thus revealing that the variational 
trial functions are not of high quality. They could be im- 
proved by adding either medium-range terms to the two- 
body Jastrow correlation or three-body Jastrow correla- 
tions. Notice that improving the importance sampling 
wave function will not affect the DMC results of Table |U 
except for the statistical error. Nevertheless, it could lead 
to better upper bounds for the L ^ energies. 

The statistical errors of the DMC energies grow 
steadily with the number of bosons. As a consequence, 
the determination of the excitation energies becomes less 
and less accurate for large N. The excitation energies El 
are obtained as the difference of two independent calcu- 
lations, one for L = (-Eoo) and the other for the desired 
value of the angular momentum (Eql). These two ener- 
gies have very close values, thus magnifying the statisti- 
cal error of their difference. As a consequence, the direct 
plot of the excitation energies will show fluctuations. In 
order to control them we have fitted all cases between 
N = 30 and N = 50 with a liquid-drop like formula, 
i.e. a third-order polynomial fit in terms of the variable 

ATl/3. 

The 3 He chemical potential, or 3 He dissociation energy, 
is defined as 



I^f = £oo( 4 H ejv ) - £oo( 4 He w 3 He) . 



(28) 



It corresponds to the energy required to eject the 3 He 
atom from the mixed system in its ground state [L = 0). 
According to the definition, hf is a positive quantity and 
its value is relevant because the states whose excitation 
energy is above it are not bound. To control the sta- 
tistical fluctuations in we have again fitted the raw 
differences with a liquid-drop formula. In Fig.n are plot- 
ted the raw DMC results for the chemical potentials and 
the excitation energies, as well as their respective fits as 
described above. For L = 4 there were too few points 
to carry out that fit, and we have only plotted the raw 
DMC results. 

The corresponding values for bound levels are dis- 
played in Table [H] for excited states with L = 1 to L = 4 
and for systems with different number of bosons. The 
total energy of the ground state (Eqo) has been already 
quoted in the last column of Table [I] The excitation en- 
ergies displayed in Table ITT1 are the raw DMC differences 
for systems with N < 20 and the results of the least 
squares fit otherwise. 

The excitation energies of radially excited levels ob- 
tained with the sum rules method described above are 
shown in Fig. [5] for L = 0,1,2. These levels are close 
to the dissociation limit, above it for N < 20 and below 
afterwards, for L = 0, 1. The radial excitation for L = 2 
is always above the dissociation limit, with the exception 
of N near 50, which signals the threshold for the binding 
of this level. One should keep in mind that the upper 
bound character of the excitation energies, as expressed 
by Eq. is not strictly satisfied in the present calcu- 
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TABLE II: Excitation energies El, in K, of He at He clusters 
for L = 1 — 4 and N — 5 — 50 in steps of 5 atoms. The 
values quoted are the result of the least squares fit described 
above for N > 20 and raw DMC results for N < 20, with the 
exception of the L = 4 column which contains the raw DMC 
results. The last column displays the 3 He dissociation limit. 



TV AE 1 


AE 2 


AE :i 


AE 4 




5 


0.48(1) 








0.53 


10 


0.37(2) 1.08(2) 






1.14 


15 


0.33(4) 0.96(5) 






1.43 


20 


0.28 


0.93 


1.64 




1.68 


25 


0.27 


0.84 


1.54 




1.85 


30 


0.25 


0.74 


1.44 


2.18(10) 


1.97 


35 


0.23 


0.65 


1.33 


2.10(12) 


2.05 


40 


0.22 


0.57 


1.22 


2.07(15) 


2.14 


45 


0.20 


0.48 


1.10 


2.01(14) 


2.21 


50 


0.18 


0.40 


0.98 


1.58(16) 


2.26 



FIG. 2: Radial excitations (in K) on top of L = to L = 2 
lowest levels. The squares are the results obtained directly 
from the DMC calculation, the continuous line is a guiding 
line for these excitations and the dashed line is the dissocia- 
tion limit. 




clearly unbound up to N = 50. 

The bound level spectra resulting from our calculations 
are collected up in Fig. 03 The bound levels are grouped 
by the value L of the angular momentum, indicated in 
the right side of the figure, with the symbol 0* signaling 
the radial excitations of the ground state. 



FIG. 3: Excitation energies of the 3 He bound excited levels 
in a mixed droplet as a function of the number of 4 He atoms. 
The levels are ordered by the value L of the angular momen- 
tum, written in the right side. Index 0* refers to the to the 
radial excitations of the ground state and corresponds in the 
figure to dashed lines. The dotted line is a guiding line joining 
the chemical potential fiF- 



2.5 r 



2- 



1.5 



0.5- 



o* 

4 



5 10 15 20 25 30 35 40 45 50 N 



10 20 30 40 N 50 




As we have mentioned in Section II, one could expect 
these spectra to look like a series of rotational bands, 
with the excitation energies roughly proportional to the 
value L(L + 1), where L is the angular momentum of the 
excited level. That is to say that the quantities 



K 



0L 



En 



L(L + 1) 



(29) 



which correspond to the rotational constants (c/. Eq.[SJ), 
should be expected to depend only on N, and not on 
the angular momentum L. This is indeed the case, as 
shown in Table lTTTl The values of the rotational constants 
smoothly decreases as N increases, as expected because 
of the dominant dependence on (l/\r F — Rb| 2 )- 



V. STRUCTURE OF 4 HEiv J HE CLUSTERS 



o t I 

10 20 30 40 N 50 



lations because the radial excitation energies have been 
computed by means of (approximate) mixed matrix ele- 
ments; therefore, some of the referred levels may not be 
bound in reality. For L > 2 the radial excitations are 



A complementary information about the nature of the 
excitations is provided by the density distributions of the 
fermion with respect to the center-of-mass of the system. 
Given that we are dealing with non-zero angular momen- 
tum states, to simplify the presentations we show them 
in Figure^] for the bound states of systems with N = 10, 
20, 30, and 40. No plot of radially excited states is given, 
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TABLE III: Rotational constants, in K, as denned in Eq. 1291 . 
for the bound levels with L = 1 — 4. 



TV 


L = 1 


2 


3 


4 


5 


0.24 








10 


0.19 


0.18 






15 


0.17 0.16 






20 


0.14 


0.16 


0.14 




25 


0.14 


0.14 


0.13 




30 


0.13 


0.12 


0.12 




35 


0.12 


0.11 


0.11 


0.11 


40 


0.11 


0.10 


0.10 


0.10 


45 


0.10 


0.08 


0.09 


0.10 


50 


0.09 


0.07 


0.08 


0.08 



because the calculation of excitation energies using sum 
rules does not allow to obtain the density distributions. 
In all the figures ps corresponds to the boson distribu- 
tion of the ground state with L = 0, as there are no 
appreciable differences between the boson distributions 
for drops in other L-excited state. This fact supports the 
model of an 3 He atom in a potential well created by the 
4 He at atoms, like in the Lekner approximation. 

Table Hvl lists the values of the root mean square radii 
of bosons and fermions with respect to the center of mass. 
Only the ground state boson radius has been displayed 
in this table since it is almost independent of the angu- 
lar momentum L. The boson radius grows monotonically 
with the number of bosons, following a rough N 1 / 3 law, 
as could be expected. The fermion radii are given for 
the states with angular momentum from L = — 4. At 
fixed L, the fermion radii increases smoothly with N, ex- 
cept for the lowest value signaling the threshold of stabil- 
ity. In that case, the fermion radius may be abnormally 
large, indicating that the system is only slightly bound. 
A flashy case is N = 30, L = 4 due to its very large 
fermion radius. This level is almost surely unbound and 
the DMC algorithm ejects the fermion far away from the 
center-of-mass of the drop; eventually the 3 He will move 
to infinity if the random walk were long enough. 



FIG. 4: The angular- averaged density distributions of the 
fermion with respect to the center-of-mass of the system 
4 Hejv 3 He for several values of the number of bosons N and 
the fermionic orbital angular momentum L. The boson den- 
sity distribution with respect to the center-of-mass appearing 
in the plot has been arbitrarily divided by 20 to appreciate 
the fermion distributions in the figures. 




10 15 



pected one. The fermion is always located at the surface 
of the boson cluster, and increasing the value of the an- 
gular momentum produces the fermion to go away from 
the boson cluster. 

To ascertain the goodness of Lekner approximation 
we have plotted in Fig. [S] the boson distributions cor- 
responding to two pure 4 He system with N and N + 1 
together with the boson distribution for N bosons plus 
one fermion. The distributions corresponding to N — 40 
are almost superimposed, thus revealing the rigidity of 
the bosonic core in front of the addition of a fermion. 
On the other hand, there are sizable differences between 
the three distributions for N = 10, revealing the effect 
of the dopant on the bosonic cluster. In other words, 
there is a weak coupling regime (Lekner approximation) 
for large N but a strong coupling regime for light drops. 



TABLE IV: Average boson (second column) and fermion 
(third to seventh columns) distances (in A) to the center- 
of-mass of the cluster. Empty entries correspond to unbound 
systems. 



N r B 


L = L 


= 1 


L = 2 


L = 3 


L = 4 


5 4.67 


6.16 


8.04 








10 5.16 


6.79 


7.87 


10.45 






15 5.49 


7.01 


8.02 


8.93 


10.74 




20 5.86 


7.50 


8.12 


9.19 


10.25 




25 6.15 


7.80 


8.42 


9.26 


10.07 




30 6.41 


8.32 


9.09 


9.67 


10.67 


17.88 


35 6.64 


8.50 


9.14 


9.87 


10.41 


11.62 


40 6.85 


8.79 


9.37 


9.90 


10.85 


11.79 


45 7.05 


9.02 


9.68 


10.18 


10.79 


11.78 


50 7.24 


9.17 


9.81 


10.47 


11.08 


12.00 



The picture which emerges from the densities plot in 
Fig. 0| and from the values of radii in Table ITVl is the ex- 



VI. CONCLUSIONS 

We have permormed DMC simulations, in conjunction 
with the sum rule method of Ref . 23] , to compute the ex- 
citation spectrum of a 3 He impurity in 4 He clusters. The 
rotational levels, namely the lowest energy levels within 
each angular momentum L subspace have been computed 
by including in the guiding function a term $i(rj?-Rfl), 
generating an eigenfunction with good angular momen- 
tum quantum numbers. On the other hand, the radial 
excitations have been estimated by computing an opti- 
mized upper bound obtained with the sum rules of order 
and 1. 

Important results have been obtained for the shell or- 
dering of the 3 He orbitals. First of all, the excitation 
spectrum contains a limited set of bound excited levels, 
whose number increases with the number N of bosons. 
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FIG. 5: Comparison of the boson distributions in A -3 , nor- 
malized to the number of particles, for the pure bosonic sys- 
tems with N and N + 1 and the mixed drop with N bosons 
and one fermion. Two cases are shown, N = 10 and N = 40. 
Labels 10 and 11 refer to pure bosonic systems, and 10+1 to 
the doped drop. In the case of N = 40 the curves are almost 
indistinguishable . 



P (A" 3 ) 
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Indicating the n-th radial excitation of the state with 
angular momentum L with notation (n + l)L, using the 
usual spectroscopic letters for the values of L, we find 
the excitation energies to follow a rotational spectrum, 
thus suggesting a shell ordering of levels Is lp Id lflg.. .. 
Starting at N w 20, the 2s radial excited state ap- 
pears as an intruder within the rotational band. Pre- 
sumably, at a larger value of TV, the 2p radial excitation 
will appear, an so on. The obtained level ordering is 
different from both the Is lp Id 2s l/2f>. . . of the three- 
dimensional harmonic oscillator and the Is 2s lp3s 2p . . . 
typical of atoms. 

This ordering of levels should be taken into account 
specially when dealing with mixed drops with a number 
of bosons much larger than the number of fermions. Note 
however that whereas the DMC algorithm is able to im- 
prove the quality of the model or importance sampling 
wave function as far as the bosonic correlations are con- 
cerned, with respect to the fermionic part it will main- 
tain the structure of the nodal surfaces. There remains 
however an important question, namely the relevance of 
these results for pure fermionic systems or for mixtures 
with a comparable number of bosons and fermions. The 
fermions are expected to play a double role: on the one 
side, fermions are creating some kind of self-consistent 
central field, analogously to the bosons, and on the other 
side they are subject to the Pauli principle effects. So, 
in a first approximation, one may assume that the level 
ordering of such systems is close to that of the system 



3 He 4 Hejv 4 +Ar 3 -i. However, one should not discard the 
analysis of other alternatives. On the basis of our results 
it becomes clear that previous microscopic calculations of 
pure fermionic drops as well as of mixed drops should be 
reconsidered, by using an improved importance sampling 
wave function based on realistic shell orderings. 
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APPENDIX A: THE SUM RULES METHOD 
APPLIED TO EXCITATIONS OF ANGULAR 
MOMENTUM L 

We have described above the use of the moment 
method to compute the radial excitation energies. The 
method may be also used to determine upper bounds to 
the excitation energies of states with angular momentum 
L, as we shall show in this Appendix. The resulting infor- 
mation will complement the one obtained directly by the 
DMC method, and is particularly relevant for the cases 
L = 1,2 where the direct calculation of the excitation 
energies is affected by a rather large relative error, being 
the difference of two large quantities, specially for large 
values of N. 

To obtain upper bounds to the excitation energy of a 
state of angular momentum L ^ it is convenient to use 
for the operator Q a form which behaves as an angular 
momentum tensor of rank L. A simple way is to consider 
the value 

Q^ L \R) = f(\r F -R B \)(x F +iy F -X B -iY B ) L . (Al) 
For the function / we have considered a power expansion 



(A2) 



with parameters C n to be determined after optimization 
of the upper bound. To fulfill the requirements which 
lead to Eq. I|13f) we should consider the Hermitian part of 
this operator. In fact, as we are interested in excitations 
of angular momentum L, irrespective of the value of the 
projection of angular momentum along some fixed axis, 
a linear combination like 



Q 



(L) _ 



(q^ + q^) /V2 



will be adequate. This combination has the advantadge 
that sum rules Mq and Mi are expressed by Eqs. I|14|) 
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TABLE V: Upper bounds to the excitation energy of the low- 
est energy state for L=l to 4 

N E x E 2 E 3 E 4 



10 0.43 1.01 1.57 2.04 1.14 

15 0.41 1.03 1.75 2.45 1.43 

20 0.36 0.92 1.59 2.27 1.68 

25 0.33 0.87 1.55 2.31 1.85 

30 0.28 0.77 1.41 2.17 1.97 

35 0.28 0.74 1.34 2.06 2.05 

40 0.25 0.69 1.29 2.06 2.14 

45 0.24 0.66 1.24 1.97 2.21 

50 0.23 0.64 1.20 1.92 2.26 



and Ijl6|l . replacing operator Q in these expressions with 

either Q( £ ) or its Hermitian conjugate . 

As in the case of radial excitations, we end up with a 
generalized eigenvalue problem, the required matrix ele- 
ments of the moment operators being given by 



M 1 
and 



2^ 



\V[r m (x-iy) L ] ■ V[r n (x + iy) L }\0) (A3) 



M° mn =(0\r m + n (x 2 +y 2 ) L \0), 



(A4) 



with r = r F - R B . Notice that < 0|Q (L) |0 >= 0. 

As the expectation values are taken with respect to the 
L = ground state, one may take the angular averages 
of the operators. After some algebra there results 



J2L-\-m-\-n 



(A5) 



and 



M x mn = ^(0\r 2L+m+n - 2 \0)[2L 2 I L ^ + {L m +Ln+mn)I L } 



with 



2fi 



I L = I [ sin 2i+1 Odd 



L\2 L 



(2L + 1) 



(A6) 



(A7) 



Particularly simple are the expressions for the upper 
bound when only the m = 0, n = matrix elements are 
retained, namely f(r) is taken as a constant, 



h 2 2(r 2L - 2 ) 1 



(AS) 



which recalls the naive the rotational model. It is worth 
mentioning the difference between this bound for angular 
excitations, Eq. I|A8(1 and the bound obtained for radial 
excitations, Eq. i|26|) . which is manifested specifically in 
the denominator. 

In Tableware given the values obtained for these up- 
per bounds, after solving the generalized eigenvalue prob- 
lem, and are quite close to the DMC excitation energies 
displayed in Table ILTI As in the case of radial excitations, 
the sum rules have been calculated by means of mixed 
matrix elements, so that they are not strictly variational. 
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